Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN WF$i style file vl.3) 



NLTE Effects in Modeling of Supernovae near Maximum 
Light 

E. Baron 1 , P. H. Hauschildt 2 , P. Nugent 1 , and D. Branch 1 

1 Dept. of Physics and Astronomy, University of Oklahoma, 440 W. Brooks, Rm 131, Norman, OK 73019-0225. 
2 Dept. of Physics and Astronomy, Arizona State University, Tempe, AZ 85287-1504- 



so 

OS 



1 February 2008 



in 



> 
oo 
oo 

o 
in 
o 

SO 

as 

6 
to 

in 
c3 



ABSTRACT 

Supernovae, with their diversity of compositions, velocities, envelope masses, and in- 
teractions are good testing grounds for probing the importance of NLTE in expanding 
atmospheres. In addition to treating H, He, Li I, O I, Ne I, Na I, and Mg II in NLTE, 
we use a very large model atom of Fe II to test the importance of NLTE processes in 
both SNe la and II. Since the total number of potential line transitions that one has to 
include is enormous (« 40 million), approximations and simplifications are required to 
treat the problem accurately and in finite computer time. With our large Fe II model 
atom (617 levels, 13,675 primary NLTE line transitions) we are able to test several 
assumptions for treating the background opacity that are needed to obtain correct UV 
line blanketing which determines the shape of near-maximum light supernova spectra. 
We find that, due to interactions within the multiplets, treating the background lines 
as pure scattering (thermalization parameter e = 0) is a poor approximation and that 
an overall mean value of e ~ 0.05 — 0.10 is a far better approximation. This is true even 
in SNe la, where the continuum absorption optical depth at 5000 A (= r st d) is << 1. 
We also demonstrate that a detailed treatment of NLTE effects is required to properly 
determine the ionization states of both abundant and trace elements. 
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1 INTRODUCTION 

Modeling the light curve and spectrum of a supernova is in 
principle a daunting task, since it requires full NLTE ra- 
diation hydrodynamics in three spatial dimensions and de- 
tailed understanding of the explosion mechanism. The dy- 
namic range and computer power for such detailed modeling 
does not currently exist. In practice one can do spherically 
symmetric, LTE, multi-group radiation hydrodynamics (e.g. 
Hoflich 1995) to determine the light curve, or one can cal- 
culate steady-state detailed model atmospheres (e.g. Baron 
et al. 1995; Nugent et al. 1995). In fact, until details of the 
explosion mechanisms are worked out and one gains the 
ability to account for complicated micro and macroscopic 
mixing, these approaches are both reasonable and accurate 
enough to extend our understanding of the nature and sys- 
tematics of supernovae, although in principle one would like 
to relax the LTE assumption in the light curve models. Even 
within this somewhat limited framework, treating all 40 mil- 
lion radiative transitions that occur in detailed NLTE is be- 
yond the capability of presently available computers, and 
there is always the caveat that the atomic data is likely to 
be inaccurate for many transitions. Nevertheless, through 



the use of a non-local approximate lambda operator we have 
been able to include (on workstation class computers) very 
detailed model atoms of iron group elements (Hauschildt 
1992, 1993; Hauschildt & Baron 1995). Our aim in this paper 
is to examine the effects of using such model atoms on the 
synthetic spectra of both Type la and Type II supernovae 
and to use the knowledge gained from numerical experiments 
to guide our necessarily approximate treatment of the the 
additional million or so radiative rates that are required to 
correctly reproduce the UV line blanketing that occurs in 
supernovae near maximum light. Actually, the task is even 
more difficult, while including about 500,000 lines saturates 
the opacity and the overall spectral shape and detailed line- 
shapes are well reproduced, which 500,000 transitions are 
the most important varies with supernova and phase and 
therefore must be dynamically determined from the full list 
of 40 million possible transitions. 

We briefly describe the Fe II model atom and the ba- 
sic features of supernova model atmospheres (details can be 
found in Baron et al. 1995; Nugent et al. 1995; Hauschildt 
& Baron 1995). Some of the description of the methods in 
Sections 2 and 3 has appeared in a companion paper on 
novae, (Hauschildt et al. 1996) although the details are pre- 
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sented in those sections for supernova models as opposed 
to nova models. The reader familiar with Hauschildt et al. 
(1996) may wish to examine the figures in those sections 
and then skip to Section 4 where we discuss results we have 
obtained for representative supernova model atmospheres. 
We conclude with a summary and discussion. 



2 METHODS AND MODELS 
2.1 Radiative transfer 

In the Lagrangian frame ("comoving frame"), the special 
relativistic equation of radiative transfer is given by (eg., 
Mihalas & Mihalas 1984) 
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Here, r is the radius, fi the cosine of the angle between a ray 
and the direction normal to the surface, v the frequency, 
I = I(r, fx, v) denotes the specific intensity at radius r and 
frequency v in the direction arccos(/i) in the Lagrangian 
frame. The matter velocity v(r) is measured in units of the 
speed of light c: /3(r) = v(r)/c and 7 is given by 7(f) = 
1/ yjl — P 2 . The sources of radiation present in the matter 
are described by the emission coefficient fj = fj(r, v), and 
X — x( r ) v ) is t ne extinction coefficient. 

We solve Eq. (1) using an operator splitting method 
with an exact non-local band-matrix approximate A- 
operator (Hauschildt 1992; Hauschildt, Storzer, & Baron 
1994). We can safely neglect the explicit time dependencies 
and partial time derivatives d/dt in the radiative transfer be- 
cause the radiation-flow timescale is much smaller than the 
timescale of the evolution of the supernova outburst (at least 
near maximum light). However, the advection and aberra- 
tion terms must be retained in order to obtain a physically 
consistent solution (Mihalas & Mihalas 1984; Peraiah 1987, 
1991a,b). This approach is also consistent with the equation 
of radiation hydrodynamics in the time independent limit 
(Mihalas & Mihalas 1984; Baron, Hauschildt, & Mezzacappa 
1996). We neither neglect the Lagrangian time derivative 
D/Dt = d/dt + vd/dr, nor assign an ad hoc value to D/Dt 
(Eastman & Pinto 1993). The latter two assumptions lead to 
physical inconsistencies with the equations of radiation hy- 
drodynamics (e.g., they do not lead to the correct equations 
for a steady-state stellar wind). Our approach is physically 
self-consistent because it includes the important advection 
and aberration terms in both the non-grey radiative transfer 
and the radiative energy equations. The only term that we 
neglect is the explicit time dependence, which is a very good 
approximation in nova and supernova atmospheres (Baron, 
Hauschildt, & Mezzacappa 1996). 

Although our method is much more complicated than 
ad hoc assumptions for D/Dt (because of the additional 
terms in the equations that must be handled that break 
the symmetry of the characteristics of the radiative transfer 
equation), its results are much more reliable than those of 



simpler methods. In addition, the solution of the correct set 
of radiative transfer and energy equations in the co-moving 
frame is no more time consuming than the solution of the 
corresponding static problem. This is because of our use of 
a non-local approximate A-operator. 



2.2 Treatment of spectral lines 

The physical conditions prevailing in nova and supernova at- 
mospheres are such that a large number of spectral lines are 
present in the line forming regions. Therefore, the simple 
Sobolev approximation cannot be used for accurate mod- 
eling including both line and continuum transfer, because 
many lines of comparable strength overlap (see Hauschildt 
et al. 1995). This means that the basic assumptions re- 
quired to derive the Sobolev approximation are not valid. 
We demonstrate this in Fig. 1 in which we plot the num- 
ber of lines that are stronger than the local b-f continuum 
and that lie within a ±2 Doppler-width wavelength inter- 
val around each wavelength point in the co-moving frame. 
This graph is for a representative Type la supernova at- 
mosphere [W7 abundances (Nomoto, Thielemann, & Yokoi 
1984), homogenized for v > 8,000 kms -1 , see e.g. Nugent 
et al. (1995)] with an effective temperature of 9,000K, and 
a micro-turbulent or statistical velocity of £ = 2kms~ 1 . In 
the UV the number of overlapping strong lines at each wave- 
length point is typically larger than 100, in some regions as 
many as 300 strong lines lie within 2 Doppler-widths. Even 
in some regions of the optical spectral range, the number 
of overlapping lines can be as high as 200 or more. Fig. 1 
shows that the situation becomes much worse for weaker 
lines (lines that are stronger than 10 -2 of the local b-f con- 
tinuum must be included in supernova atmosphere models) . 
Now, nearly 1000 lines (all of comparable strength) overlap 
in the UV and around 300 lines in the optical can overlap 
at each wavelength point. This shows decisively that the 
Sobolev approximation cannot be used in modeling super- 
nova atmospheres. In light of this observation, our approach 
solves the full radiative transfer equation for all lines. In ad- 
dition, we wish to model stellar atmospheres (both static 
cool stars, and hot stars with winds) and in those cases the 
Sobolev approximation does not apply. 



3 THE FE II MODEL ATOM 

We have taken the data needed to construct our model atom 
from the compilation of Kurucz (1994), whose long term 
project to provide accurate atomic data for modeling stellar 
atmospheres is an invaluable service to the scientific com- 
munity. For our current model atom, we have kept terms up 
to the 2 H term, which corresponds to the first 29 terms of 
Fe II. Within these terms we include all observed levels that 
have observed b-b transitions with log gf > —3.0 as NLTE 
levels (where g is the statistical weight of the lower level and 
/ is the oscillator strength of the transition) . This leads to a 
model atom with 617 levels and 13675 "primary" transitions 
which we treat in detailed NLTE. We solve the complete b-f 
and b-b radiative transfer and rate equations for all these 
levels and include all radiative rates of the primary lines. 
In addition, we treat the opacity and emissivity for the re- 
maining nearly 1.2 million "secondary" b-b transitions in 
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Figure 1. The lower line shows the number of overlapping lines within two Dopplcr widths of each wavelength point when lines are 
included that are at least as strong as the bound-free continuum. The upper line shows the same when weak lines that are within 0.01 
of the bound-free continuum are included. 



NLTE, if one level of a secondary transition is included in 
our detailed model. 

Using this procedure to select our model atom, we ob- 
tain 13675 primary Fe II NLTE lines between the 617 Fe II 
levels included in NLTE. For Fe II lines with A > 3500A, 
we use 5 to 11 wavelength points within their profiles. In 
extensive test calculations, we have found that Fe II lines 
with A < 3500A do not require these additional wavelength 
points in their profiles due to the fine wavelength grid used in 
the model calculations at these wavelengths (Hauschildt & 
Baron 1995). This procedure typically leads to about 30,000 
wavelength points for the model iteration and the synthetic 
spectrum calculations. Figure 2 illustrates the adequacy of 
this procedure. We compare the output spectra of two mod- 
els, which were iterated independently. One has the line by 
line opacity sampling below 3500A and the other has full line 
profiles for all NLTE lines. Clearly this wavelength grid is 
sufficiently dense for the calculations. This is mostly due 
to the properties of the Fe II ion, in particular that its 
lines are concentrated mostly in the wavelength range be- 
low 3500A. Other model atoms, e.g., Ti I, require additional 
wavelength points for practically every line (Hauschildt et 
al, in preparation). For all primary lines the radiative rates 
and the "approximate rate operators" (see Hauschildt 1993) 
are computed and included in the iteration process. A de- 
tailed description of the Fe II line treatment is presented in 
Hauschildt et al. (1996). 



4 RESULTS 

In this section we compare and contrast the effects of NLTE 
Fe II in SNe II, where iron is only a trace element with 
those in SNe la where iron is a major constituent. We will 
discuss the transition class SNe Ib/c in future work. Table 1 
describes the basic model parameters for the spectra that 
we use to illustrate the NLTE effects. 

4.1 SNe II 

In SNe II, iron is only a trace element in the outer ejecta, 
and thus, one might expect that its effects on spectrum for- 
mation would be small. This is not the case, since iron line 
blanketing still plays the dominant role in the formation of 
the spectrum. Figures 3 and 4 display calculations with Fe II 
treated in LTE and NLTE compared to observed spectra of 
SN 1987A 13 and 31 days after explosion. The varying role 
of Fe II NLTE is nicely displayed in these figures. While we 
have attempted to fit the observed spectrum, neither the 
LTE nor NLTE models match the width of the all of the ob- 
served features, since composition inhomogeneities as well 
as more complicated density structures (including clump- 
ing) are known to be important in SN 1987A. Nevertheless, 
the NLTE Fe II models overall, particularly in the blue, pro- 
vide a better fit to the observed spectra than do the LTE 
Fe II models. 

In Fig 3 the Ha and H/3 absorption is too narrow in 



both models and the absorption feature at 5800A, which is 
likely due to He I A5876 and the Na D doublet, is poorly 
fit. However, below A ~ 4200A, the NLTE model more 
closely follows the observed spectrum. In the later spectrum 
(Fig. 4), the overall fit of the NLTE model is again supe- 
rior to that of the LTE spectrum, even though the Balmer 
lines are now quite complicated probably due to composition 
inhomogeneities . 

4.2 SNe la 

In SNe la, where iron is a major constituent of the atmo- 
sphere, one expects that Fe II NLTE effects would play a 
dominant role. Figure 5a compares the results of two repre- 
sentative SNe la model calculations where Fe II is treated 
in LTE and NLTE. In both cases the models have been iter- 
ated to radiative equilibrium. In order to better display the 
differences between the two calculations, Fig. 5b displays 
the differences between the two spectra, where a positive 
value implies more flux in the NLTE model. The differences 
are significant, particularly since one goal of spectral analy- 
sis through synthetic spectra is to determine abundances in 
supernovae. 

Up to this point when we have referred to LTE treat- 
ment, we have really meant transitions, for which the rates 
are not included in the solution of the rate equations. How- 
ever, the source function is not treated in complete LTE in 
that we write the source function in the line as: 

S*a = eB x + (l-e)J x , 

where e is the thermalization parameter (which we take to 
be a constant in the range 0.05 < e < 0.1) and the other 
symbols have their usual meaning. In a pure LTE treatment 
e = l, and in a "pure scattering" treatment e = 0. Figure 6 
compares the results of a full NLTE treatment, a pure LTE 
treatment, and a pure scattering treatment of the lines, with 
the structure held fixed. The results are quite interesting. 
The pure LTE case reproduces the overall spectral shape of 
the NLTE case rather well, while the pure scattering case 
does not. In fact the overall luminosity differs in the NLTE 
case and the pure scattering case by more than a factor of 2. 
The lineshapes, are however, better reproduced in the pure 
scattering case. This result is not surprising, since in many 
lines an analytic estimate of e such as given by the equiva- 
lent two level atom formulation would predict a very small 
value of e < 0.01. On the other hand, the total shape of 
the spectrum depends on the interaction between lines and 
continua as well as the interaction among overlapping lines 
(i.e. blanketing, blending, and collisions in the multiplets). 
The dotted line shows a case where only 594 b-b transitions 
have been included and illustrates that the "pure scatter- 
ing" approximation does not correctly pass flux from blue 
to red. Thus, even though it reproduces the lineshapes in 
the optical, and since it has the opacity to correctly reduce 
the flux due to line blending, the flux is simply scattered 
and not passed to longer wavelengths as it should be. In 
the NLTE case, the multitude of collisions within multiplets 
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Figure 2. The solid line shows the result of a calculation of an SN la model with every line profile function is included with 5 points 
per line for a total of 83269 wavelength points and the dashed line shows the results of sampling for A < 3500 A, with a total of 38465 
wavelength points. 



Table 1. Parameters of Presented Models 



Model 


^"rriodcl 


Ro 






VQ 


^Ni 


Abund 


N 


SN 87A_day_13 


5400 


6.1 x 10 14 


cm 


5500 


km s _1 


0.03 


0.33 Z e 


6 


SN 87A.day.13 


5100 


4.5 x 10 14 


cm 


1700 


km s _1 


0.01 


0.4 Z Q 


6 


SN la 


9000 


9.7 x 10 14 


cm 


7500 


km s — 1 


0.05 


W7 





The basic model parameters for the spectra presented. The model atmospheres are characterized by the following parameters (see 
Hauschildt et al. 1992; Baron et al. 1995, for details): (ii) the model temperature T mo< j e i, which is defined by means of the luminosity, 
L, and the reference radius, Ro, (T mo( jei = (L / '4-77 Ro 2 a) 1 ' ' 4 where a is Stefan's constant); the reference radius Ro, which is the radius 
where the continuum optical depth in extinction at 5000 A (r st( j) is unity; the expansion velocity, ^o, at the reference radius; the mass 
fraction of 56 Ni, X^i, which is the source of non-thermal gamma-rays; the elemental abundances which are for the SN 1987A models 
solar (Anders & Grevesse 1989) with mctalicity scaled by the factor given in the table and the SN la model are W7 (Nomoto, 
Thielemann, & Yokoi 1984) compositions, homogenized for v > 8,000 kms -1 . the density structure parameter, v e , 

(p(r) oc cxp(-v/v e )), or N (p(r) oc (v/v )~ N ) 



(see also Hoflich et al. 1996) creates a pseudo-thermal pool 
of photons which are in equilibrium among themselves and 
have a more nearly isotropic momentum distribution. 

4.3 NLTE effects on the Fe II ionization 

In Fig. 7 we display the ground state departure coefficients, 
b\, for two supernova models. Since supernovae produce sig- 
nificant amounts of radioactive 56 Ni it is important to in- 
clude the effects of non-thermal electron (produced by 7-ray 
deposition) collisions in the NLTE rate equations. In the 
models presented here the 7-ray deposition is treated as lo- 
cal due to a nickel mass fraction Xi-h, described in Table 1. 
The effects of non-thermal collisional ionization by primary 
electrons produced by collisions with gamma rays due to 
the decay of 56 Ni and 56 Co are modeled using the continuous 
slowing down approximation (Garvey & Green 1976; Swartz 
1991; Baron et al. 1996). We have also included ionizations 
due to secondary electrons (Swartz 1991; Baron et al. 1996), 
but most of their energy is thermalized and therefore has 
only a small effect on the level populations (Meyerott 1980). 
The collisional cross sections are taken from the work of Lotz 
(1967a,b, 1968a,b,c). 

Figures 7a,b show the results for the model SN la with 
and without non-thermal effects, whereas Fig. 7c, d shows 
the results for the model SN 87A_day_13, with and without 
non-thermal effects. We have marked the bi for Fe II with 
connected plus-signs for clarity. In all models, the departures 
from LTE are significant. In most regions of the supernova 
atmospheres, the Fe II 61 are smaller than unity. This usu- 
ally indicates an overionization of Fe II relative to LTE, but 
we will see below that this is not always the case. The de- 
parture coefficients of other species show the same behavior 
as Fe II, which is very different from the results we have ob- 
tained in our nova atmosphere modeling (Hauschildt et al. 
1996). The sharp drop of the departure coefficients of He I 
in the model SN 87A_day_13 is caused by the ionization of 
helium by non-thermal electrons in this model as can be seen 
by comparing the models with and without non-thermal ef- 
fects. Most effected by the non-thermal effects are He I, Fe II, 
and Mg II. Ca II appears to be hardly affected. Comparing 
the SN la models with and without non-thermal effects, it 



is evident that the non-thermal effects play an important 
role in determining the level populations. Again Ca II, ap- 
pears to be the least affected by non-thermal effects. In the 
SN II model O I is not strongly affected by non-thermal col- 
lisions although that is not the case in the SN la model. In 
the SN la model, NLTE effects are quite small when non- 
thermal ionization is neglected. This may be due to the fact 
that since SN la atmospheres consist entirely of metals, with 
lower ionization energies than hydrogen and helium, that 
thermal electron collisions act to restore LTE. Once all im- 
portant non-thermal electron collisons are included NLTE 
effects become important. 

If Fig. 8 we display the ionization structure of iron 
throughout the atmosphere for model SN 87A_day_13. In 
panels (a and b) we show the NLTE results (with and with- 
out non-thermal electrons) whereas panels (c and d) show 
the results for LTE (but otherwise the same model struc- 
ture). First note the number of iron ionization stages that 
are present in the atmospheres (in a number of interleaved 
"Stromgren Spheres" ) . This is a more pronounced and well- 
known feature of nova atmospheres which is caused by the 
large temperature gradients inside the expanding shell. We 
find similar effects for all elements that are included in the 
model calculations. 

The effects of NLTE on the Fe II ionization balance be- 
come clear by comparing panels (a and c) and (b and d) in 
Fig. 8. In the LTE case, Fe II is the dominant ionization stage 
of iron in the outer atmosphere, In the NLTE-Fe II case, 
however, Fe II recombines into Fe I throughout the outer at- 
mosphere, although the ground state departure coefficients 
of Fe II are less than unity (the latter fact usually indi- 
cates a overionization relative to LTE) . This is true whether 
or not non-thermal effects are included, although the iron 
remains ionized to lower optical depths when non-thermal 
effects are included. The reason for this non-intuitive ioniza- 
tion behavior becomes clear by inspecting the graphs shown 
in Fig. 9. In the NLTE case (Fig. 9a), hydrogen remains 
ionized throughout the outer atmosphere, and the electron 
densities therefore remain high. However, in LTE, hydrogen 
recombines and the electron densities are much smaller than 
in the NLTE model. This completely changes the behavior 
of the ionization equilibrium for all species and the higher 
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Figure 3. The observed spectrum of SN 1987A (taken at CTIO) on March 8, 1987 (thin solid line) is compared to two model calculations 
with reference radii and velocities of 6.1 X 10 14 cm and 5500 kms -1 and T mo< j e i = 5400 K, X^; = 0.03, and Z = 0.33 Z@. The density 
profile is a powerlaw with N = 6 (model SN 87A_day_13 in Table 1). The thick solid line includes Fe II in NLTE, whereas the thin 
dot-dashed line displays a calculation where Fe II is treated in LTE. 



Figure 4. The observed spectrum of SN 1987A (taken at CTIO) on March 26, 1987 (thin solid line) is compared to two model calculations 
with reference radii and velocities of 4.5 X 10 14 cm and 1700 kms -1 and T mo d c i = 5100 K, Xni = 0.01, and Z = OAZq. The density 
profile is a powerlaw with AT = 6 (model SN 87A_day_31 in Table 1). The thick solid line includes Fe II in NLTE, whereas the thin 
dot-dashed line displays a calculation where Fe II is treated in LTE. 



electron pressures combined with the different recombina- 
tion rates in the NLTE model force the recombination of 
Fe II to neutral iron. This effect illustrates that NLTE mod- 
eling is both extremely important and must be done in great 
detail in order to obtain realistic results. In particular, this 
shows that to fit an observed SN II using a model in LTE, 
would force the choice of a high temperature in order to re- 
produce the Balmer lines. With this high temperature, in 
LTE the ionization of Fe could never be reproduced. The 
correlation of the recombination of iron with the ionization 
of hydrogen is evident upon comparing Fig. 9a and Fig. 9c. 
In panel (a) which includes non-thermal electron collisions, 
the iron is neutral deeper into the atmosphere, whereas in 
panel (c) without non-thermal effects the iron is only recom- 
bining in the outer regions just where the electron and H II 
fractions are rising. 

In Figure 10 we show that the changes in the iron ion- 
ization structure in the SN I model are very small. This 
is expected since the ground state departure coefficients for 
Fe II are close to unity in this model. Again, with the neglect 
of non-thermal effects, the high electron fractions in SN la 
models drive them towards LTE. This is very different from 
the behavior we find in both the SN II model and our nova 
model atmospheres (Hauschildt et al. 1996). The NLTE ef- 
fects of Fe II in the SN I model will be apparent only in the 
line source functions and the spectra, as discussed below. 

In Fig. 11 we show the departure coefficients for all Fe II 
levels for the two models (with and without non-thermal ef- 
fects) in an overview graph. This figure demonstrates that 
the NLTE effects for Fe II are important and differ for the 
individual levels and for different model temperatures, (for a 
direct comparison of effect of varying the model temperature 
in SNe II, see Hauschildt & Baron 1995). For each T mo( jei, the 
departure coefficients are closer to unity (their LTE value) 
for the higher lying levels, which are more strongly coupled 
to the continuum than the lower lying levels. The rather 
large variation of the bi within the Fe II model atom suggests 
that the Fe II lines will show significant NLTE effects them- 
selves, in addition to the effects introduced by the changes 
in the ionization structure. Comparing the models with and 
without non-thermal effects we see that the models with- 
out non-thermal ionization have departure coefficients larger 
than unity at moderate optical depths, however when non- 
thermal ionization is considered the departure coefficients 
tend to stay below one. 



4.4 NLTE effects on the formation of Fe II lines 

The NLTE effects not only change the ionization balance of 
iron but also change the relative level populations of Fe II. 



We demonstrate this effect for two important Fe II multi- 
plets. In Figs. 12a,b we show the ratio of the line source 
functions Sl to the local Planck function for the UV1 mul- 
tiplet of Fe II. These are UV lines (around 2600A) from 
the ground a 6 D term to the z 6 D° term. In the model SN la 
which includes non-thermal effects (Fig. 12a, left panel) the 
UV1 multiplet is nearly in LTE up to the outer atmosphere. 
Only below r std = 1CF 3 do the line source functions dif- 
fer, by about 20%, from the Planck function. The situation 
is different for the model SN 87A_day_13 (Fig. 12a, right 
panel). Here, the line source functions are larger than the 
Planck function by nearly a factor of 2. In both models, 
the (collisional) coupling between the levels of each term 
is strong enough to nearly establish the same line source 
function to Planck function, Sl/B, ratio for all lines within 
the multiplet, although the electron densities are relatively 
low. Only in the SN II model do the levels within a term fi- 
nally decouple for r stt j < 10~ 4 . When non-thermal processes 
are neglected, the same general trend is observed although 
now the SN II model has a UV1 line source function that is 
smaller than the local Planck function inside, but becomes 
greater than the local Planck function far out in the atmo- 
sphere. 

The subordinate Fe Il-multiplet 42 (a 6 S to z 6 P°) shows 
a somewhat different behavior in Fig. 12a. The line source 
functions for the 3 primary NLTE lines are slightly smaller 
than unity for 10~ 4 < r st( j < 2 and larger than unity for 
T st d < 10~ 4 for the model SN la (left panel). The model 
SN 87A_day_13 (right panel) shows nearly the opposite be- 
havior. The NLTE effects are slightly larger for the cooler 
model, as seen earlier for the UV1 multiplet. This shows that 
the NLTE effects for Fe II are not just over-ionization but 
that the "internal" NLTE effects are also significant. This 
multiplet is less affected by non-thermal ionization than is 
the UV1 multiplet. 

The Fe II NLTE effects we have found in SN I and II 
atmospheres are significantly smaller than in nova atmo- 
spheres (Hauschildt et al. 1996). This is due to the much 
higher velocity field in SN atmospheres, which tends to 
transfer photons from the regions neighboring the Fe II lines 
into them, thus coupling the frequencies stronger together 
than in the case of novae. This is also true for the regions 
around the ionization edges of Fe II, which limits the effec- 
tive number of Fe II ionizing photons. It must be noted here 
that our departure coefficients, radiation field and source 
functions are all computed in the co-moving frame and, 
therefore, simplified arguments based on Eulerian frame ex- 
perience can and will lead to wrong conclusions. In the co- 
moving frame, the velocity does not lead to a frequency shift 
of the photons (by its construction) but it leads instead to 
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terms which are better described as "drag" terms, slowing 
the response of the radiation field to changes in the opac- 
ity as a function of wavelength when compared to the static 
case. This "stiffer" radiation field limits the magnitude of 
the NLTE effects on the line formation, in fact this is more 
effective at higher velocities. Therefore, stellar winds and no- 
vae show comparatively larger NLTE effects than the very 
fast expanding atmospheres of SNe. 



5 CONCLUSIONS 

We have demonstrated that large model atoms are required 
in order to correctly transport energy and push the flux from 
blue to red in SNe atmospheres. Weak lines must be included 
in order to correctly reproduce line blanketing and blending. 
In spite of the fact that an equivalent 2-level atom approach 
would predict very small values of the thermalization pa- 
rameter e for many transitions, collisions in the multiplets 
and continuum effects redistribute the energy and momen- 
tum such that an overall higher value of e leads to a more 
accurate description of the radiation transport. A calcula- 
tion that treats weak lines as pure scattering misses far too 
much opacity and gives results similar to those that treat 
only continua. Accurate modeling of the UV also requires 
large model atoms and the inclusion of weak lines as back- 
ground opacity. As better atomic data become available the 
accuracy of the calculations should improve. NLTE effects 
cause very non-intuitive ionization effects to occur which 
show that these effects must be included in detail. 
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Figure 5. Panel (a):The thick solid line corresponds to a SN la model where Fe II is treated in NLTE, and the thin solid line corresponds 
to a model where Fe II is treated in LTE. Both models are iterated to radiative equilibrium and have the parameters: W7 (Nomoto, 
Thielemann, & Yokoi 1984) compositions, homogenized for v > 8,000 kms^ 1 , reference radii and velocities of 9.7 X 10 14 cm and 
7500 kms -1 and T mot j e i = 9000 K, X^i = 0.05 (model SN la in Table 1). The density profile is a exponential in velocity with scale 
v e = 900 kms -1 . The thick solid line includes Fe II in NLTE, whereas the thin dot-dashed line displays a calculation where Fe II is 
treated in LTE. Panel (b):The relative difference between the spectra shown in Panel (a), where a positive value corresponds to more 
flux in the NLTE model. 

Figure 6. The results of a model calculation where the lines are treated in full NLTE (thick solid), in pure LTE (thin solid), in pure 
scattering (short-dashed), and where only the 594 b-b transitions of the model atoms Mg II, Ca II, O I, Na I, Ne I, H I, He I, and He II 
are included (long-dashed line). 

Figure 7. The run of the ground state departure coefficients bi for the NLTE species as functions of T st( j (the b-f optical depth at 
5000A). Panels (a,b) show the results obtained for the model SN la, with and without non-thermal ionization, respectively and panels 
(c,d) the model SN 87A_day_13, with and without non-thermal ionization respectively. 



Figure 8. Ionization structure of iron as function of T st( j (the b-f optical depth at 5000A) for the model SN 87A_day_13. Panel (a) shows 
the results obtained with Fe II treated in NLTE, with non-thermal ionization, panel (b), NLTE without non-thermal ionization, whereas 
panel (c) shows the results for LTE Fe II, with the structure of the non-thermally ionized model, and panel (d) shows LTE Fe II, with 
the structure of the model neglecting non-thermal ionization. The NLTE undcr-ionization causes Fe II to recombine to Fe I in the outer 
layers of the SN atmosphere. Note the number of "Stromgren Spheres" of various ionization stages of iron in the inner parts of the SN 
atmosphere. 



Figure 9. The concentrations of the most abundance species as function of T s td (the b-f optical depth at 5000A) for the model 
SN 87A_day_13. Panel (a) shows the results obtained with in NLTE, including non-thermal ionization, panel (b), NLTE without non- 
thermal ionization, whereas panels (c, d) show the results for LTE, respectively. The NLTE under-ionization causes Fe II to recombine 
to Fe I in the outer layers of the SN atmosphere. Note the number of "Stromgren Spheres" of various ionization stages of iron in the 
inner parts of the SN atmosphere. Non-thermal processes prevent this recombination from occuring deeper in the atmosphere. 



Figure 10. Ionization structure of iron as function of r st( j (the b-f optical depth at 5000A) for the model SN la. Panels (a, b) show 
the results obtained with Fe II treated in NLTE, with and without non-thermal ionization, respectively, whereas panels (c, d) show the 
results for LTE Fe II, with and without non-thermal ionization respectively. In this model, the changes of the ionization of iron are very 
small and visible only in the outermost parts of the atmosphere. Note the number of "Stromgren Spheres" of various ionization stages 
of iron in the inner parts of the SN atmosphere. 



Figure 11. The departure coefficients for all NLTE levels of Fe II as functions of the standard optical depth r st( j for two models: SN la, 
with and without non-thermal effects (panels a and b) and SN 87A_day_13, with and without non-thermal effects (panels c and d). The 
departure coefficients are closer to unity for the higher lying levels. 



Figure 12. The ratio of the line source function to the local Planck function for the UV1 multiplet (a 6 D to z 6 D°) and for multiplct 
42 (a 6 S to z 6 P°) as functions of the standard optical depth. The left side displays the model SN la and the right side the model 
SN 87A_day_13. Panel a includes non-thermal ionization and panel b does not. 



